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Abstract 

We show that Plebanski's second heavenly equation, when written as a first- 
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Therefore by Magri's theorem it is a completely integrable system. Thus it 
is an example of a completely integrable system in four dimensions. 
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1 Introduction 

The Einstein field equations that govern self-dual gravitational fields reduce 
to a single scalar- valued equation. This is either the complex Monge- Ampere 
equation that Calabi pQ had shown to govern Ricci-flat Kahler metrics and 
which Plebanski |2| has called the first heavenly equation, or 

u u u xx - u tx 2 + u xz + u ty = , (1) 

which is Plebanski's second heavenly equation. In this paper we shall con- 
sider the Hamiltonian structure of the second heavenly equation. We shall 
show that it can be formulated as a Hamiltonian system in three or more 



inequivalent ways. Therefore, by the theorem of Magri j3j it is a completely 
integrable system. The only known example of a completely integrable sys- 
tem in four dimensions was anti-self-dual Yang-Mills fields. Here we show 
that the second heavenly equation which governs anti-self-dual gravitational 
fields is another example of a completely integrable system in four dimen- 
sions. 

The real Monge-Ampere equation admits bi-Hamiltonian structure j3] 
and in the case of the complex Monge-Ampere equation we have two sym- 
plectic structures jS] except in the crucial case of two complex dimensions. 
It is natural to expect that the second heavenly equation (JIJ may admit bi- 
Hamiltonian structure. This expectation is further supported by the general 
feeling that self-dual gravitational fields should be integrable systems. 

Earlier Dunajski and Mason [S| mentioned briefly that Plebanski's sec- 
ond heavenly equation may admit bi-Hamiltonian structure. For this purpose 
they rewrote (0) as a nonlocal first order equation. They used a scalar re- 
cursion operator appropriate to the original one-component second heavenly 
equation. They did not explicitly construct Hamiltonian operators determin- 
ing the structure of Poisson brackets, nor the family of commuting Hamil- 
tonians. In their later paper [7j Dunajski and Mason suggested another 
recursion operator for the one-component equation ((!} determined by two 
different, though compatible, recursion relations but they did not develop 
the Hamiltonian formalism there. 

We shall consider a different, a two-component representation of the sec- 
ond heavenly equation as a first order system of two local equations and show 
that Plebanski's second heavenly equation indeed admits Magri, i.e. multi- 
Hamiltonian, structure. Earlier we had constructed [5] the scalar recursion 
operator for the second heavenly equation and now we shall cast it into 2x2 
matrix form which naturally joins the two recursion relations into one matrix 
relation. 

In section El we introduce a first-order two-component form of the second 
heavenly equation . In section |3] we present the first Hamiltonian struc- 
ture of this system of equations. We start with a degenerate Lagrangian and 
construct its Dirac bracket jH] to find the Hamiltonian operator. In section 0] 
we invert the Hamiltonian operator and obtain the corresponding symplec- 
tic 2-form. In section El we construct explicitly a matrix integral-differential 
recursion operator in the two-component form which incorporates naturally 
both recursion relations. This operator and the operator determining sym- 
metries form a Lax pair for the two-component system. In section El we give 
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explicitly the second and the third Hamiltonian structures which show the 
way for obtaining multi-Hamiltonian representation of our two-component 
system. In section |7] some simple symmetries of our equations are shown to 
be generated by certain integrals of motion via previously constructed Hamil- 
tonian operators. We find the complete Lie algebra of point symmetries of 
the two-component system and for all variational symmetries we construct 
the corresponding integrals of motion. In section[H]we give examples of higher 
flows obtained with the aid of the Hermitian conjugate of the recursion op- 
erator. Those flows are nonlocal symmetries generated by local integrals. 



2 First order form of the second heavenly 
equation 

The second heavenly equation is a second order partial differential equation. 
In order to discuss its Hamiltonian structure we shall single out an inde- 
pendent variable, t, in (0) to play the role of "time" and express the second 
heavenly equation as a pair of first order nonlinear evolution equations. Thus 
we introduce an auxiliary variable q whereby ([!} assumes the form 

u t = q , q t = — (q x ~%~ u xz ) = Q (2) 

of a first order system. For the sake of brevity we shall henceforth refer to 
(J2J) as the second heavenly system. It is worth noting that this split of (0) 
into the system (J2J) is not unique, here we are using the most straight-forward 
choice. The choice of the independent variable t as "time" is arbitrary and 
in no way represents any sort of a physical time variable. Also we shall 
henceforth assume u xx ^ 0. This is not an essential restriction but rather a 
statement of non-triviality. Now the vector field 

defines the flow. In the discussion of the Hamiltonian structure of this system 
we shall use matrix notation with 

u l i = 1,2 u 1 = u, u 2 = q 

running over the dependent variables. 
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The equations of motion (J2J) are to be cast into the form of Hamilton's 
equations in two different ways according to the recursion relation of Magri 



u\ = X(u l ) = J l k 6 k H! = Jf 5 k H (4) 

where S k denotes variational derivative of the Hamiltonian functional with 
respect to u k . 

3 First Hamiltonian structure 

There is a systematic way to derive the first Hamiltonian structure of (J2J). 
It was used to obtain Hamiltonian structures for the real Monge-Ampere 
equation We shall now apply it to the second heavenly system. 

We start with an action principle for (J2J given by the Lagrangian density 

111 

C = qu t u xx + -u t u y - -q 2 u xx + -u x u z (5) 

which is degenerate. We need to apply Dirac's theory of constraints |Hj in 
order to arrive at its Hamiltonian formulation. We define canonical momenta 

dL 




which satisfy canonical Poisson brackets 

[*i{t),Av)] = 8t6(t-T,) 

where £, rj are generic names for independent variables each one of which 
stands for the collection of our original independent variables x, y, z. In other 
words, 8(£ — r)) = 5{x — x') S(y — y') 6(z — z') for £ = {x, y, z}, rj = {x', y', z'}. 
We note that the momenta cannot be inverted for the velocities because the 
Lagrangian (j£J) is degenerate. Therefore following Dirac [5] we impose them 
as constraints 

<Pu = 7T U - (qu xx + ~ Uy) (6) 

and calculate the Poisson bracket of the constraints 

K ik = [<pi(x, y, z), <p k (x\ y\ z')\ (7) 
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organizing them into the form of a matrix. We find 



[<f) u (x,y,z),(f) u (x',y',z')] 

[<f) u (x,y,z),(f) q (x\y',z')] 
[$ q (x,y,z),<f> q (x',y',z?)] 



- q(x) 5 x i x i(x - x) S(y- y')5(z - z') 
+ q(x) 5 xx (x - x') 5{y - y')5{z - z') 
-5(x-x')5 y (y-y')5(z-z') (8) 
-u xx 8{x- x') 8{y- y') 6(z - z') 




and see that they do not vanish modulo the constraints. Thus the constraints 
are second class in the terminology of Dirac. We need to invert the matrix 
of Poisson brackets of the constraints in order to arrive at the Dirac bracket. 
That is why we have not simplified (jS)) using properties of the Dirac delta- 
function jTUJ. We shall do so in sequel, (fTo")) . but it is better to keep them in 
raw form for purposes of inversion because the definition of the inverse 



j K ik (Z,a)J k i(a,r ] )da = 5i5(Z-r 1 ) 



(9) 



results in a set of differential equations for the entries of J fej . 

Given any two smooth functions of the canonical variables A, B the Dirac 
bracket is defined by [0] 



m),B(v)] D = [A{t),B( V )] 



(10) 



where J lk is the inverse of the matrix of Poisson brackets of the constraints 
defined by ©. 

The first Hamiltonian operator for second heavenly system is given by 
the Dirac bracket 



/ 



jik 
J 



1 



J_ \ 



V 



^D x + D x -% 

XX XX 



(11) 



and it can be directly verified that 

1 



^i = 2 ^ Uxx ~ 2 UxUz 



;i2) 



is the conserved Hamiltonian density for the flow (J2J). In other words, the 
total t-derivative of Tii along the flow (J2J) is a three-dimensional total diver- 
gence in "space" variables 

dHi _ dP 1 dQ 1 dRt 
dt dx dy dz 

and hence the total t-derivative of the functional H\ = f^TCidxdydz van- 
ishes by the divergence theorem, provided the components of the current 
Pi, Qi, Ri approach zero fast enough as x, y, z go to ±00. 

The Hamiltonian functional obtained from the density (|12|) is the one 
that yields the equations of motion (J2J) by the action of the operator (fTTj). 

The proof of the Jacobi identities for the Hamiltonian operator (fTTj) is 
straight-forward but rather lengthy. A shorter proof is obtained by inverting 
(fTTj) to arrive at the symplectic structure of the second heavenly system. 



4 Symplectic structure 

The statement of the symplectic structure of the equations of motion (j2J) 
consists of 

i x u = dH (13) 

which is obtained by the contraction of the closed symplectic 2-form uj with 
the vector field X (jSJ) defining the flow. The symplectic 2-form is obtained 
by integrating the density 

u = -du 1 A Kijdu j (14) 

where K is the inverse of Jq given by 1)1 ip. But we already know it from the 
Poisson brackets of the constraints (JHJ). Thus the inverse of the Hamiltonian 
operator (fTTj) is given by 

K= [ Uxx ^ 



uj = q x du A du x — u xx du A dq — -du A du y (16) 



which is a local operator. Then we find the symplectic 2-form from (fPf)) 

1 

— i 
2 

which, up to a divergence, (fTo^) can be directly verified to be a closed 2-form. 
The closure of the symplectic 2-form (fTB*|) is equivalent to the satisfaction of 
the Jacobi identities for the Hamiltonian operator (jllj). 
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5 Recursion operator 



Recently a recursion operator for second heavenly equation was obtained 
[7j, [S]. We shall use it to construct new Hamiltonian operators satisfying 
Magri's recursion relation (J3J). But before we can do so we need to express 
the recursion operator in the two- component form appropriate to the system 

We start with the equation determining the symmetries of the second 
heavenly system. We introduce two components for symmetry characteristics 

UT = iP (17) 



q T = ip 

of the system (J2J). From the Frechet derivative of the flow we find 

A-( Dt 2 ~ l ) (18) 

\ "ra "'xx "'xx "ra " / 

and the equation determining the symmetries of the second heavenly system 
is given by 

A($) = 0. (19) 

We note that the combination of the first and second determining equations 
f!19|) with the latter multiplied by an overall factor of u xx , coincides with the 
determining equation for symmetries of the original second heavenly equation 
(JTJ). It was noted in jS] that this determining equation has the divergence 
form 

(u xx i[> - q x ip x + ip y )t + (q t f x ~ Qx^P + ¥z)x = (20) 

rewritten in two-component form. This implies the local existence of the 
potential variable (p such that 

fit = qt¥x - q x ip + 

<fix = -{Uxxlp ~ q X fx + <fy) (21) 

which, as is proven in j^j, satisfies the same determining equation for sym- 
metries of (JTJ and therefore is a "partner symmetry" for (p [TT]. In the 
two-component form we define the second component of this new symmetry 
similar to the definition of ip as ijj — (p^. Then the two-component vector 




satisfies the determining equation for symmetries in the form (|19|) and hence 
is a symmetry characteristic of the system (J2J) provided the vector (|17|) is 
also a symmetry characteristic. Thus (|2*Tj) become the recursion relation for 
symmetries in the two-component form 

$ = K(<f>) (22) 

with the recursion operator given by 

-p _ ( D x l {QxD x - D y ) -D~ l u xx \ . , 

K ~ \ QD X + D z -q x ) {26} 

where D~ x is the inverse of D x . See |T2] for the definition and properties of 
this operator, in particular, 

D ^=l(L-D mdt (24) 

and the integrals are taken in the principal value sense. The commutator of 
the recursion operator ()23|) and the operator determining symmetries (fTHj) 
has the form 

/ D- 1 (q t -Q) xx -(q t -Q) x D~ 1 (u t -q) X i 

[n,A]= [ {^ t -Q) xx + (D y -^)(q t -Q) x }D x (q t -Q) x 

(25) 

and as a consequence, the operators 1Z and A commute 

[K, A}=0 (26) 

by virtue of the second heavenly system (j2J). Moreover, TZ and A form a Lax 
pair for the second heavenly system. 



6 Second and third Hamiltonian structures 



The second Hamiltonian operator J x is obtained by applying the recursion 
operator (J2*3|) to the first Hamiltonian operator J\ = 7ZJ . We find 



( D~ l 

X 



Jl 



U, 



V'XX 



/ qx_ y) _1_ _|_ _1_ y) _£x_ _|_ _i_ n + £j _1_ 

2 V ^J^XX ^ ^^XX ^^XX ^ ^^XX '^'XX ^ ^* '^'XQ 
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(27) 



which is manifestly skew. The proof of the Jacobi identity is again straight- 
forward and lengthy. 

The Hamiltonian operators (fTTJl and (|27j) form a Poisson pencil, that is, 
every linear combination a Jo + bJ\ of these two Hamiltonian operators with 
constant coefficients a and b satisfies the Jacobi identity. This can be verified 
using the functional multi- vectors criterion of Olver ^3] • 

The vanishing of the total time derivative of the Hamiltonian functional 

= J^°oo 'Hodxdydz with the Hamiltonian density 

H = (x + Ci)qu xx , (28) 

where c\ is an arbitrary constant, follows by the divergence theorem from the 
fact that &Ho/dt, calculated along the flow (|2*|. is a three-dimensional total 
divergence in x, y, z and therefore Ho is an integral of the motion along the 
flow (J2J), provided the components of the current tend to zero fast enough as 
x, y, z approach ±00. 

The Hamiltonian function 7i satisfies the recursion relation (@J) 

u\ = JpSkH! = J\ k 5 k H (29) 

which shows that the second heavenly equation (P) in the two-component 
form (J2J) is a hi- Hamiltonian system. 

In constructing the second Hamiltonian operator we have used the fact 
[Sj that given one Hamiltonian operator Jo and the recursion operator 1Z 

Jn = n U J (30) 

is also a Hamiltonian operator. In the case of (|27j) we have n — 1. Now if 
we act with the recursion operator (J23j) on the second Hamiltonian operator 
Ji, or use (|3Up for n = 2 we can generate a new Hamiltonian operator J 2 = 
IZJi = J\ Jq 1 Ji where we have used the fact that by construction 1Z = J\ Jq 1 . 
The explicit expression for J2 is 

1 -D- l D y D-' -( D -i Dz „%+^.\ \ 



qx(q y + u X2 )D x \- + -l-D x q x (q y + u xz ) 

XX XX 

\ u xx w xx 

- i qxDz ^x~x + ltx- Dzq * 



(31) 
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which is again manifestly skew and we have checked the Jacobi identity for 
J 2 and its compatibility with J\ and J . 

It can be verified that the Hamiltonian functional with the density 

H-i = u - (y + c 2 ) qu xx , (32) 

with C2 another arbitrary constant, is an integral of the flow (J2J). Its im- 
portance follows from the fact that it also satisfies the Magri's relation (Jlj) 
with Ji replaced by the third Hamiltonian operator J2, so that we obtain 
a tri- Hamiltonian representation of the second heavenly equation (JTJ in the 
two- component form (J2J) 

u] = 4%H X = Jl%H = Jf 5 k H_ v (33) 

This construction can be continued by the repeated application of the recur- 
sion operator (|23|) to Hamiltonian operators. Thus we have the Magri, or 
multi- Hamiltonian representation of second heavenly system (J2J). 

7 Symmetries and integrals of motion 

Hamiltonian operators provide a natural link between commuting symmetries 
in evolutionary form and conserved quantities, integrals of motion, in 
involution with respect to Poisson brackets. Our original two-component 
system (J2J) is also a member of the infinite hierarchy of commuting symmetries 
and has the form fl33J) where H u H and H-i, defined bv (fT2j) . and 
respectively, are integrals of motion, while the whole right-hand side of (}3*3*j) 
is the symmetry characteristic of time translations. If we replace Jo by J\ 
and Ji by J 2 on the right-hand side of an d the time t by the group 
parameter r we obtain 

< = Jf = Jf 4/Jo = «*, (34) 

that is, the symmetry of ^-translations generated either by the integral Hi 
relative to the second Poisson structure J x , or by the integral H relative to 
the third Poisson structure J 2- We also note that 
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= -( u xU y - u 2 x q x ). (37) 



where C3 is another arbitrary constant and furthermore 

= -(T) (36) 

where c<i is an arbitrary constant coming from H_i. 

Now we consider the functional H 1 determined in a similar way to Hi by 
the density 

1 

1 

and check that the total time derivative of Ti 1 along the flow (J2J) is a total 
divergence and hence the functional H l is an integral so that it can serve as 
a Hamiltonian for some flow commuting with Acting on the column of 
its variational derivatives by Jo and Ji we obtain two such flows 

< = J$5 k H l =u x 

< = jj^if 1 = 

with symmetry characteristics of x and ^-translations respectively. 

Now let us present the results of a systematic search for point symmetries 
and the corresponding integrals of motion for the second heavenly system. 
The complete Lie algebra of point symmetries of the original one-component 
heavenly equation (0) was given in [Hj but in the two-component representa- 
tion the results look a little different. The basis generators of one-parameter 
subgroups of the complete Lie group of point symmetries for the second 
heavenly system (J2J) have the form 

Xi = -2zd t + txd u + xd q , X H = td t + zd z + ud u , W f = f(y, z)d u 

X IH = td t + xd x + 3ud u + 2qd q , Z b = b{y)d z - b'(y)xd t - b"{y)— d u 

o 

/ tx 2 x 2 \ 
Y a = ad y + a'(xd x - td t - zd z + qd q ) + a" [ xzd t — d u - — d q ) 



3 

m x z 




+ a'" — d u , V d = d z (y,z)(td u + d q ) - d y (y,z)xd u (38) 

U c = Cyd t + c z d x - c yz x(td u + d q ) + c yy — d u + c 2 

where a(y), b(y), c(y, z), d(y, z) and f(y, z) are arbitrary functions, primes de- 
note ordinary derivatives of functions of one variable and we used the short- 
hand notation d t = d/dt and so on. Since some of the generators contain ar- 
bitrary functions, the total symmetry group is an infinite Lie (pseudo)group. 
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In the table of commutators of the generators (|38|) the commutator 
stands at the intersection of ith row and jth column and we have used the 
shorthand notation / = zf z — /, / = af y — a'zf z and B(y) = J b(y)dy. Here 
= c y d z — c z d y is the Jacobian. 





x I 


X n 


Xjii 


Ya 


z b 


Uc 


v d 


w f 


X I 








Xi 





u 2B 


—V2C+C 


-w 2zdz 





X H 














-z b 


U c 


v d 


Wf 


Xrn 


-Xj 
















-2V d 


-3W f 


Y a 














Z(ab)' 


U~c 


v d 


Wf 


z b 


-u 2B 


z h 





—Z(ab)t 





u bCz 


V bdz 


w bfz 


u c 


V26+C 




u c 


-u~ c 


-u bCz 





Wd(c,d) 
S(m,z) 





v d 


w 2zdz 


-y d 


2V d 




-v bdz 


— Wd{c,d) 








w, 





-Wf 


3W f 


-Wf 


-w hf . 












Table 1: Commutators of point symmetries of the second heavenly system. 

We complete this table with the commutation relations [Y a , Y g ] = Y ag t^ ga i, 
[Z b , Z h ] = 0, [Uc, U s ) = V aM , [V d , V e ) = and [W f , W r ) = where g = g{y), 

d(y,z) 

h = h(y), s = s(y, z), e = e(y, z) and r = r(y, z) are arbitrary functions. 

We are interested in the integrals of motion generating all these point 
symmetries. The relation between symmetries and integrals is given by the 
Hamiltonian form of Noether's theorem 




Jo 



5 U H 
S q H 



(39) 



where H = J_™ Tidxdydz is integral of the motion along the flow (J2J), with the 
conserved density 7i, which generates the symmetry with the two-component 
characteristic JT^] Vu,Vq- We choose here Poisson structure determined by our 
first Hamiltonian operator J since we know its inverse K given by (|15|) which 
is used in the inverse Noether theorem 



6 U H 




(40) 



determining the integral H corresponding to the known symmetry fj u ,rj q . 

We proceed by using the formula (|40|) for reconstructing conserved den- 
sities corresponding to all variational point symmetries. For the symmetry 
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X I we use its characteristic fjj = (xt + 2zq,x + 2zQ) T , with T denoting 
transposition, in the formula (J4(J|) to obtain the conserved density 



Hi = (txq + zq )u xx + - u x - zu z u x . 



x 



(41) 



For the symmetries Xjj and Xjjj corresponding integrals in (j4*Uj) do not exist 
and hence they are not variational symmetries jT3|. A characteristic of the 
symmetry Y a has the form 

/ ^a"'x 3 z — ^a"tx 2 — (a"xz — a't)q — a'xu x — au y + a'zu z \ 
^ a \ a'q — ^a"x 2 — (a"xz — a't)Q — a'xq x — aq y + a'zq z J 

and from ()40|) a conserved density corresponding to it is given by 



Ha 



q 



+ u 



(ta' — xza") + q 



—x za — -tx a + a \zu z — xu x ) — au y 



Un 



\ t x 2 a!" - -x 3 za IV 
2 6 



1 „ t 1 



+ -a xu x yzu 
+ -a'(zu y u z - tu x u z - xu x u y ) - ^au 2 y 



-xu r 



(42) 



where a(y) is an arbitrary smooth function. A characteristic of the symmetry 
Z b is 



Vb 



— \x 3 b" + xb'q — bu z 



xb'Q — bq z J 
and from ()4()|) we obtain the corresponding conserved density 

„3 \ 



Jib 



X -b'q 2 -{^b" + bu z \q 



u xx + ^x 3 b"'u - ^{b'xu x + buy)u z (43) 



where b(y) is an arbitrary smooth function. A characteristic of the symmetry 
U c is 



Vc 



tc zz xCy Z CyQ c z q x 



and the corresponding conserved density calculated by (|4*U|) has the form 



(a - c z u x )q - c 



y 



1 

— ( 

2 



- OyU + -a t u x + -(c y u z - c z u y )u x (44) 
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where c(y, z) is an arbitrary smooth function and we have used the notation 



( \ 2 ^2 

(T[t) X, y, Z) —t C ZZ tXCy Z ~\- —X Cyy. 

A characteristic of the symmetry Vd is given by rjd = {td z —xd y , d z ) T and the 
corresponding conserved density is 

ft<z = (td z - xd y )qu xx - (td yz - xd yy )u + ^d z u 2 x (45) 

where d(y, z) is an arbitrary smooth function. Finally, the symmetry Wf has 
the characteristic fjf = (f, 0) T and the corresponding conserved density is 

ft/ = fqUxx - fyU (46) 

where f(y, z) is an arbitrary smooth function. 

By a lengthy calculation one may check that the time derivatives of all 
these densities Ji along the flow (J2J) are total divergences providing an in- 
dependent check that the corresponding functionals H are indeed integrals 
of motion subject to suitable boundary conditions. Note that by replacing 
q by Ut we obtain integrals of motion for the original form ((TJ of the second 
heavenly equation. 

Finally, we list some simple obvious symmetries, which are particular 
cases of the symmetries Y a , Z b , U c and Vd, and integrals of motion corre- 
sponding to them. Among the latter we find all the Hamiltonian functions 
7ii,Ho,H_i and TC 1 (up to inessential arbitrary constants) which were used 
so far. Translational symmetries in each of the four coordinates and integrals 
of motion generating them are 

ft* = ft c =j/ = -^(q 2 Uxx + u x u z ) = -Hi 

(47) 

qu x u xx + -u x u y ^j H x = -- (u x u y - u 2 x q x ^j = -H 1 

H y = H a =i = - (qu y u xx + ^u 2 y 
W = Hb=x = - (qu z u xx + ^u y u z 



x t 


= d t = Uy 


x x 


= d x = U z 


n x 


= H c=z = 


Xy 


= dy = Y 1 


x z 


= d z = z l 
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The dilatational symmetry appears as the particular case of Y a 

X d = Y y = xd x + yd y - td t - zd z + qd q (48) 
H d = H a =y = ^(q 2 u xx - u x u z ) + (zu z - xu x - yu y ) (qu xx + -u^j . 

Boosts in z- and x-directions appear as the special cases of Z b and U c 

X B * =Z y = yd z - xd t (49) 

H Bz = H b=y = [^q 2 - yqu z ^j u xx - ^{xu x + yu y )u z 

X B * = U yz -X I = zd t - yd x 

U Bx = - Qg 2 - yqu x ^j u xx + -(yu y + zu z )u x 

with the latter symmetry combined with Xj. Note that the symmetries 
X/, Xji and X in are not particular cases of Y a , Z b , U c , Vd and Wf. Shifts in 
u and q appear as special cases of Vd and Wf 

x st = td u + d q = v z n st = n d=z = t qUxx + x -u 2 x 

X Sx = xd u = V d =- y H Sx = H d= - y = xqu xx = H (50) 

X Sy = -yd u = Wf=^ y H Sy = Hf=- y = -yqu xx + u = H-i. 

Note that the second heavenly equation itself has divergence form 

(u x q x - Uy) t = (u x q t + u z ) x (51) 

and therefore h = u x q x — u y , that is 

h = Hf=-i = -qu xx (52) 

is a conserved density and can therefore serve as a Hamiltonian. To find 
the Hamiltonian flows, we apply the Hamiltonian operators Jo and J\ to 
the vector of variational derivatives of the Hamiltonian functional C = 
I-oo h d x dy dz corresponding to h 

(Sg)-(i) < 53 » 
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with the results 



and 



and hence 



5 U C 
5 9 C 




5 fc C = 0, n > 1. (55) 

The last formula shows that C is the Casimir functional relative to the Pois- 
son structure operators J n with n > 1, that is, these operators have a non- 
trivial kernel and therefore are non-invertible on the whole phase space, as 
opposed to Jq, which implies non-invertibility of the recursion operator as 
well. Hence the phase space equipped with any of the Poisson structures 
Jn, n > 1 is a Poisson manifold but not a single symplectic leaf. 



8 Higher flows 

We know from the work of Fuchssteiner and Fokas [H] (see also ^H] and 
references therein) that if a recursion operator has the form TZ = JiJ ' 1 , 
where Jo and J\ are compatible Hamiltonian operators, then it is hereditary 
(Nijenhuis), i.e. it generates an Abelian symmetry algebra out of commuting 
symmetries. Moreover, Hermitian conjugate hereditary recursion operator 
acting on the vector of variational derivatives of an integral yields the vector 
of variational derivatives of another integral. Therefore we calculate TV to 
find 

| (D x q x -D y )D-' -D X Q-D Z ^ 



and act by it on the vector of variational derivatives of H\ with the result 

( u z q xx - q z u xx + 2q x u xz - u 





(57) 



T<2 = q u z u xx - ~u u yz (58) 
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is a new integral. This also can be checked straightforwardly by computing 
the total time derivative of Ti 2 along the flow (J2J) and showing that it is a 
total divergence. 

By definition TZ = JiJq 1 we have TV = Jo Ji and therefore (joTjl can be 
rewritten as 

J*6 k H 2 = jf 4#i = <. (59) 

A nontrivial result is obtained by applying J\ to the vector of variational 
derivatives of H 2 which yields 

( u r\ =J [ b u H 2 \f D- 1 D z {u x q x -u y )-u x q z ^ 
[q T J 1 ^ S q H 2 J ~ y Q Uxz + Uzz _ qxqz 

where the right-hand side is a nonlocal symmetry characteristic. In this case 
the local integral H 2 generates a nonlocal symmetry. 

We obtain similar results using H 1 with the density ()37|) instead of Hi. 
The action of TV on the vector of variational derivatives of H 1 yields again 
a vector of variational derivatives 

-rjf ( \ I u y q xx — q y u xx + 2q x u xy — u yy \ / 5 U H 2 \ . . 

K \ 6 q Hi ) ~ { u y u xx ) ~ { 5 q H* ) {bl} 

where 

H 2 = qu y u xx + -u 2 y (62) 

is a new conserved density. Acting by Jo on the variational derivatives 
of H 2 will not yield new results since V) = Jq Ji and hence J^5kH 2 = 
—Ji SkH 1 = u y . Using Ji instead of Jo yields 

j ( 6 ^H 2 \ ( D- l D y (u x q x - u y ) - u x q y \ 
\ s i H J \ Uyz- q x q y + UxyQ j 

which is a nonlocal symmetry characteristic generated by a local integral H 2 . 

In a similar way we can construct higher integrals and corresponding 
higher flows by applying the conjugate recursion operator TV to the varia- 
tional derivatives of all the integrals constructed in section 
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9 Conclusion 



We have cast the second heavenly equation into the form of a two-component 
local nonlinear evolutionary system in order to discover its Hamiltonian struc- 
ture. We started by presenting its first Hamiltonian structure. Then we cast 
the scalar recursion operator for the second heavenly equation into matrix 
form. Thus we were able to obtain explicitly the second and third Hamilto- 
nian structures for the second heavenly system. The recursion operator and 
the operator determining symmetries form a Lax pair for the two-component 
system. By Magri's theorem the multi-Hamiltonian structure makes the 
second heavenly system a completely integrable system in four dimensions. 
Apart from anti-self-dual Yang-Mills this is first truly completely integrable 
system in four dimensions. 

Following Magri, we define the (complete) integrability of Hamiltonian 
equations as the existence of an infinite number of conservation laws and 
symmetries related to them by the Poisson structure or, equivalently, the 
existence of a recursion operator for symmetries. To justify the name 'inte- 
grable' we should show how this property will lead to obtaining solutions, 
that is to the integration of the system. In our papers [TTJ |S| we suggested 
the method of partner symmetries for obtaining solutions of the complex 
Monge- Ampere equation and the second heavenly equation of Plebanski. The 
existence of partner symmetries was closely related to the existence of the 
Lax pair of Mason and Newman ff)J but there is no known version of 
the inverse scattering method in four dimensions which could utilize this Lax 
pair in a customary way. A particular choice of partner symmetries provided 
differential constraints such that the equation under investigation together 
with constraints possessed an infinite Lie group of contact symmetries which 
resulted in a linearizing transformation and exact solutions. Thus we can 
consider the method of partner symmetries as an alternative use of Lax pairs 
for linearization of the problem as opposed to inverse scattering method. Our 
new method for integrating four- dimensional heavenly equations is still in its 
early stages of development and it is not the subject of our study in this 
paper. We plan to return to this matter in the future. 

We have presented the Lie algebra of all point symmetries and the inte- 
grals of motion which generate all variational point symmetries of the second 
heavenly system. We gave examples of some nonlocal symmetries (higher 
flows) generated by local Hamiltonian functions through these Hamiltonian 
structures. We are in the process of applying a similar approach to the 
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complex Monge- Ampere equation. 
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